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Abstract 

The gap equation in the 1 5'o-channel is solved for a nuclear slab with the separable form 
of the Paris potential. The gap equation is considered in the model space in terms of the 
effective pairing interaction which is found in the complementary subspace. The absolute 
value of the gap A turned out to be very sensitive to the cutoff K max in the momentum 
space in the equation for the effective interaction. It is necessary to take K max =16Q 4- 
180 fm _1 to guarantee 1% accuracy for A. The gap equation itself is solved directly, 
without any additional approximations. The solution reveals the surface enhancement of 
the gap A which was earlier found with an approximate consideration. A strong surface- 
volume interplay was found also implying a kind of the proximity effect. The diagonal 
matrix elements of A turned out to be rather close to the empirical values for heavy 
atomic nuclei. 
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1 Introduction 



The microscopic theory of nuclear pairing is yet far from being complete. Some progress was 
achieved only recently. In a series of papers (see pQ and Refs. therein), the problem was 
attacked within a two-step approach in which the full Hilbert space of the problem S is split 
into the model subspace So and the complementary one, S'. The gap equation should be 
solved in the model space in terms of the effective pairing interaction (EPI) V^ s which obeys 
the Bethe-Goldstone-type equation in the complementary subspace. The pairing effects could 
be neglected in the latter provided the model space is sufficiently wide. The simplest nuclear 
systems with one-dimensional (ID) inhomogeneity were considered, i.e. semi-infinite nuclear 
matter and the nuclear slab. These simplifications permitted to find the EPI directly without 
additional approximations for the separable representation E] of the Paris NN-potential [I] . 
To save the computing time, it was suggested also an approximate method of finding the EPI, 
the so-called local-potential approximation (LPA), which has very high accuracy PQ. 

At the same time, the gap equation in P was solved only approximately, with the help of a 
generalization to non-homogeneous systems of the method by V.A. Khodel, V.V. Khodel and 
J.W. Clark (KKC) jSJ. Here we report the result of the direct solution of the gap equation for 
the nuclear slab system. 

In recent quite relevant paper by F. Barranco et al. jB], the direct solution of the gap 
equation for a finite nucleus ( 120 Sn) was carried out for the realistic Argonne v 14 NN-potential. 
It yielded the gap A which is about a half of the experimental value. The lack of the gap was 
attributed in jB] to the contribution of the low-lying surface vibrations. Comparison of the 
results of our calculation with those of Ref. [H] will show to what extent the predictions of the 
microscopic calculation of A in finite systems depend on the form of the realistic NN-potential 
and other details. In particular, it is well known [7j that the gap value in infinite matter is 
rather sensitive to the value of the effective mass m*. It is illustrated with Fig. 1 for the case 
of the Paris NN-potential. Qualitatively, such a strong dependence of A could be explained by 
the simplest BCS formula for A in the case of the momentum independent force: 



where g is the dimensionless effective pairing interaction strength, and the energy coefficient Eq 
depends on the cutoff in the integral over momentum in the BCS gap equation. 

In [6. the coordinate dependent effective mass m*(r) was used of the Sly4 version jHj of the 
Skyrme force. Its value changes from m* = 0.7m inside a nucleus to m* = m outside. We prefer 
to deal with the prescription of m* = m which is closer to the "experimental" nuclear value. 
The latter was found in Ref. 9j from adjusting the single-particle spectra of magic nuclei, in 
agreement with the older analysis by G. E. Brown et al. [TUj . 

2 Outline of the formalism 

Following to pQ , we use the many-body theory form of the gap equation fTJ 




(1) 




(2) 
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which explicitly takes into consideration the particle-particle propagator in the superfluid sys- 
tem, 

A'(E,e) = G(j+e)G'(j-e). (3) 

Here G and G s stand for the single-particle Green functions in normal and superfluid systems 
correspondingly. Obvious space and spin variables in Eqs. (J2J) and (j3J) are omitted. As usual, 
the product in Eq. (j2J) means the coordinate integration and spin summation. The total two- 
particle energy is E = 2fi, where /i is the chemical potential under consideration. 

Within the Bethe-Brueckner approach, the irreducible particle-particle interaction block V 
is approximated by the free energy-independent iViV-potential. In this case, the gap A is also 
energy-independent and the r.h.s. of Eq. (j2J) can be integrated over e' explicitly yielding the 
well-known expression of the gap in terms of the anomalous density matrix k: 

/c(ri, r 2 ) = J A S (E; n, r 2 , r 3 , r 4 )A(r 3 , r 4 )dr 3 dr 4 . (4) 

Here we use the same as above notation A S (E) for the integral of Eq. Q over e. In a 
symbolic notation, the gap equation reads: 

A = -Vk. (5) 

When the solution of the gap equation in a finite system is carried out in a direct way, the 
form (jSJ is preferable as far as the anomalous density matrix is expressed directly in terms of 
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the Bogolyubov t>-functions as follows: 

fe(ri, r 2 ) = Mi(ri)fi(r 2 ), (6) 

i 

where the summation is carried out over the complete set of the Bogolyubov functions with 
energy eigenvalues Ei > 0. 

The form (J2J) of the gap equation is more convenient for the two-step method of solution 
described in the Introduction. After the splitting of the complete Hilbert space S = So + S' 
described above, the two-particle propagator is represented as the sum A s = + A'. Here we 
already neglected the superffuid effects in the S"-subspace and omitted the upperscript "s" in 
the second term. 

The gap equation (J2J) can be rewritten in the model subspace: 

A = -KMA, (7) 

where V^ s obeys the following equation: 

V? s = V - VA'V&. (8) 

The model space S is defined by a separation energy E in such a way that it involves 
all the two-particle states (A, A') with the single-particle energies e x ,e x > < E . In this case, 
the complementary subspace S' involves such two-particle states for which one of the energies 
e\, E\i or both of them are greater than E Q . 

The gap equation ((7j) in the model space can also be written in the form 

A = -V>o, (9) 

where 

k = A S A. (10) 

However, a simple expression similar to Eq.fjUJ) doesn't exist for the anomalous density 
matrix k in the model space. Indeed, let us substitute the usual pole expansions for Green 
functions G, G s JT] for the integral over e of the expression similar to Eq. © for the model space 
propagator Aq. Let us then expand the u-, f-functions and the gap A in the basis of the single- 
particle functions 4>\ (r) which makes the normal Green function G diagonal, G\\>(e)=G\(e)5\\i. 
One gets: 

u*(r) (11) 

A 

«i(r) = X;« < VA(r), (12) 

A 

A(n,r 2 ) = £A AA ,0 A (n)<Mr 2 ). (13) 

A, A' 

After simple operations, one finds 



,0! / n Xl v^v^ ( (l-n Xl )^u^ 




K [T 1 ,T 2 )= > U + 



Ei - £ Xl + \l Ei+E Xl - yU J 

x A AlA3 Al (n)0 A2 (r 2 ). (14) 
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Here e\ are the single-particle energies without pairing and Ei, the ones with pairing taken into 
account. The upperscript "(0)" in the sum means that all the states A belong to the model 
space: E\ < E . 

Let us now write down explicitly the Bogolyubov equations in the "A-representation" : 

A' 

£ A AV uf - (e x - (jl) = Ei vf , (15) 

A' 

With the use of Eqs. ()15j) . the sum in Eq. (JT%j) can be simplified to the following form: 

(o) 



k 



r 1; r 2 )= Yl (nxM'^ + il-nxJufv^Mn^x^). (16) 



i.X-t ,A 2 



It can be easily seen that in the limit of Eq — > oo, when the set (A) is complete, the expression 
(fTBj) coincides with Eq. (JSJ). 

Following to [I], we use the separable 3x3 form [21 E| of the Paris potential [I] in the 1 Sq 



channel, 



V(k,k') = J2\j9 i (k 2 )9 J (k' 2 ), (17) 



where the form factors gi{k 2 ), (i = 1, 2, 3) are given in PQ. 

In this paper, we deal with a slab of nuclear matter embedded in a potential well U(x), 
where x is the direction of the lD-inhomogeneity. In this case, the momentum in the (y, z)- 
plane (or s-plane) is the integral of motion. The single-particle wave functions can be written 

as 

A (r) = ^ s y n (x), (18) 
and the Bogolyubov functions have the form: 

u n (v) = e l ^ s u l (k ± ,x) (19) 

v i (r)=e ik ± e v i {k±,x), (20) 

In this case, it is natural to use the mixed representation, i.e., the coordinate representation 
for the x-direction and the momentum one for the s-plane. In this representation, the separable 
NN-potential leads to a similar form of which, in the notation of pQ, reads: 

V c p s (k 2 ± , k'l Xl , x 2 , x 3 , x 4 ; = £ Aij(X, X'- E) 9i (kl, x)g 3 (k'l x'). (21) 

Here the center-of-mass and relative coordinates in the x-direction are introduced (X = (xi+X2)/2, 
x = X1—X2, etc.), and gi(k\,x) stand for the inverse Fourier transform of the form factors 
gi(k\ + k 2 ) in the x-direction. Their explicit form is given in pQ. 
The coefficients Ay obey the set of integral equations: 

A ij (X,X';E) = \ ij 8(X-X')+ 

+ J2 X ^J dX 1 B lm (X,X 1 ;E)A mj (X 1 ,X';E), (22) 

lm 



where B\ m are given by the convolution integrals of the propagator (—A') with two form factors. 
Their explicit form is as follows: 



B lm {X,X'-E)=Y! J 



d 2 k ± (1 - N k±>n - N k±y ) 
(2tt) 2 E-e n - e n , - k 2 Jm 

xg l nn ,(ki,X)g™(kl,X% (23) 



stvn>(k±,X) = J dx gi (kl,x)y n (X+x/2)y n ,(X-x/2). (24) 

In Eq. (|23|) Nf- ±jn =(0, 1) are the occupation numbers in the system without pairing. The 
symbolic summation over (nn') denotes the summation over the discrete spectrum for negative 
energies and the integration over dpdp' / (2ir) 2 , for positive ones. All the momentum integrals in 
()23|) are limited with a cutoff, p,p', k ± < K max . The prime in the sum of Eq. shows that 
the summation is carried out over (nn') which are not included in the model space of Eq. (JJJ) 
for A. In the above short notation, we have (AA') G S'. Therefore, the energy denominator 
in Eq. f!23[) always exceeds the difference (Eq — p) ^> A what permits to neglect the pairing 
contribution to this propagator. 
Analogously to Eq. (J2IJ, we have 

A(k ± ;x u x 2 ) = J2MX)g l (k 2 ±1 x) 1 (25) 

i 

and 

K (k±; x 1 ,x 2 ) = Yj ^o( x )9i(kl, x). (26) 

i 

The gap equation (JJJ) is reduced to the set of three one-dimensional integral equations, 

Ai(X) = " E / A ^'( X > X '-> <{X') dX'. (27) 

3 

For the ID-geometry, the Bogolyubov functions can be expanded in the basis of y„-functions: 
Ui(k±,x) = ^2ui(k±)y n (x) , Vi(kj_, x) = ^2v™(k±) y n (x) . (28) 

n n 

Then, for a fixed value of k±, the Bogolyubov equations have the form: 

(e n - p) u?(k ± ) +J2 A nn'(k ± ) vf(k ± ) = E iU ?{k ± ) , 
n' 

Y,^nn>(k ± ) u?(k ± ) - (e n - p) v?(k ± ) = E iV ?(k ± ) , (29) 



where 



Then we have 



A nn ,(k ± ) = Y, j M*) 9nn^X) dX. (30) 

i 

= 1^1 E Ko n '(ki_)gi n ,(ki,X) , (31) 

J \ Z7V ) nn' 



where 

< n '(k ± ) = £ (N k± ,n u?(*j_) vf(k ± ) + (1 - iV fc± , n ) u»'(* ± ) v?{k ± )) • (32) 

i 

The chemical potential /i in is, as usual, determined by fixing the average particle 
number N. In the system under consideration, one should fix the number of particles per a 
unit surface of the slab: 

°- = / E «(^)) 2 = const - ( 33 ) 

The above set of equations is complete and can be solved by an iteration procedure. 



3 Calculation results and discussion 

In calculations, we used, just as in [I], the one- dimensional Saxon- Woods potential U(x) sym- 
metrical with respect to the point x = 0, 

U<yX) = 1 + exp((x - L) /d) + exp(-(x + L) /d) ' (34) 

with typical for finite nuclei values of the potential well depth Uq = 50 MeV and diffuseness 
parameter of d = 0.65 fm. For the main part of calculations the thickness parameter was chosen 
as L = 8 fm to imitate heavy nuclei. The dependence of the gap A on the value of L is also 
examined. 

It should be noted, that a simplification of the numerical procedure in the slab system under 
consideration can be made using the parity conservation which follows from the symmetry of 
the Hamiltonian under the axis reflection x — * —x. As a result, the eigenfunctions y n can 
be separated into even, and odd, y~, functions. Then the two-particle propagators in 
Eqs. ©j© split into the sum 

A = A + + A- (35) 

of the even and odd components. The first one, A + , originates from the terms of the sum in 
Eq. containing states (A, A') with same parity, and the second one, A", from those with 
opposite parity. So long as the NN-potential V does conserve the parity, the propagators A + 
and A~ do not mix in Eq. (|5|). Therefore the separation similar to Eq. (|H5|) holds true for the 
correlation part of the EPI, 6V^ S —V^ e — Vq. 

As far as we deal with the singlet 1 So-pairing, the gap A is the even coordinate function. 
Therefore only even components of the 2-particle propagators Aq, A' and of the correlation part 
5V? S enter Eqs. © and ©• 

Let us now discuss several computational details. To simplify account of the continuum, we 
used the discrete spectrum representation method by putting the infinite wall at the distance 
R > 30 fm. Independence of results on the value of R was checked. The iteration procedure of 
solving the gap equation was stopped when the maximum difference between values of Ai(X) 
found in the consecutive iterations became less than 10~ 5 MeV. The typical number of the 
effective iterations was 40. At each iteration, the chemical potential \x entering the Bogolyubov 
equations was found anew from Eq. (jHEJ)- 
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By solving the gap equation, one finds a set of the components Aj(Jf), 2, 3 of the gap 
function and the ones, Ki(X), of the anomalous density matrix. To present the results in a 
more transparent way we introduce, following to P, the "Fermi averaged" gap 

A F (X) = J2&i(X)9i(4(X)), (36) 

i 

and density matrix, 

K F (X) = Y / K l (X)g l (kUX)). (37) 

i 

Here the local Fermi momentum is defined as follows: k F {X)=2m{ix— U(X)) if 2m(fi—U(X)) > 
0, and kp(X)=0 otherwise. For the Fermi averaged EPI, we have 

V e F ff (X) = J2 hAXMkl{X))gM{X)l (38) 

ij 

where the zero moments of the EPI components Aij(Xi, X 2 ) are defined in a standard form: 

Aij(X) = J dtA ij (X-t/2,X+t/2). (39) 

We found that the model space with the maximum energy E =20 MeV is sufficiently wide 
for convergence of the calculation procedure. It is illustrated in Fig. 2. 
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Figure 2: The Fermi averaged gap A-p(X) for different E . 

It is seen that the function A F (X) changes noticeably under variation from Eq=0 to 
E =10 MeV and only a little, from E =10 MeV to E =20 MeV. Finally, it doesn't practi- 
cally change under variation from Eq=20 MeV to £'o=30 MeV. The Fermi averaged EPI at 
different E is displayed in Fig. 3. 
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Figure 3: The Fermi averaged EPI Vj s (X) for different E { 



It reveals a strong attraction outside the slab which changes rapidly in the surface region 
to rather weak attraction inside the slab. The asymptotic value of Vf s (X — > oo) decreases 
with growing of E . As E — > oo, the function V^ S (X) must tend to the analogous quantity Vq 
calculated for the free NN- interaction. The asymptotic value of is Vq(X — > oo) = — 720 MeV. 

The situation with convergence of the integrals of Eq. ()23j) for the 2-particle propagator 
turned out to be more dramatic. A very slow convergence of this integral was noted in 
where the problem of microscopic calculation of the EPI was analyzed, but the gap equation 
itself was not solved. Such a big value of the cutoff momentum as K max =60 fm -1 was chosen to 
guarantee a reasonable accuracy of the EPI. The main reason of this phenomenon is the hard- 
core character of the Paris potential. Now we found that the gap A is much more sensitive to 
the cutoff than the EPI, and such a huge value as K max =160 -j- 180 fm" 1 is necessary to obtain 
A with 1% accuracy. It can be seen from Fig. 4 for the Fermi averaged gap Ap(X) calculated 
for different values of K max . The corresponding results for the EPI are displayed in Fig. 5. 

To emphasize the difference between the curves, we displayed in Fig. 6 the regions of small 
and large X separately, in a enlarged scale. It is seen that the variation of the EPI with growing 
of K max from 60 fm -1 to 180 fm -1 looks quite small, but it leads to a noticeable enhancement 
of A-p(X). In fact, 1% variation of EPI results in 5%, or even more, change of A. The reason 
of such a high sensitivity of A to a small variation of the EPI could be seen in the exponent 
dependence of A on the inverse effective interaction strength 1/g in the BCS formula (0). In 
this approximation, all depends on the value of g. Experience of solving the gap equation 
with Paris potential in infinite nuclear matter shows that in the region of the maximum of A 
(k,F — 0.8 fm -1 ), where the effective g is big, the convergence is not so slow. In this case, 
K max =20 fm -1 is enough to reach 1% accuracy in A. However, let us consider the density 
describing the inner region of the slab under consideration, kp=lA2 fm -1 , which corresponds 
to a small A (small g) situation. In this case, the convergence is very slow, and for 1% accuracy 
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Figure 4: The Fermi averaged gap Ap(X) for different cut-off momenta 
Kmax- The model space corresponds to E = 20 MeV. 
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Figure 5: The Fermi averaged EPI Vf s (X) for different cut-off momenta 
K m ax- The model space corresponds to E = 20 MeV. 

in A one needs in K max ~ 200 fin . Looking in Fig. 5, one could conclude that the surface 
region with a big EPI should dominate in the gap equation ()27j1 . In this case, the sensitivity 
of A to small variation of the EPI should not be so high. However, the comparison of Fig. 4 
and Fig. 5 shows that the inner region with small EPI (and high sensitivity of A to g) should 
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Figure 6: Details of Fig. 5 in an enlarged scale. 

contribute to Eq. (J2Tj) significantly to explain this strong influence of a variation of K max on 
the value of A. 

To check this point, it is instructive to analyze the Fermi averaged anomalous density matrix 
kf(X). It is displayed in Fig. 7 together with the EPI V^ S (X) and the product of Vf s (X) Kp(X). 
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Figure 7: The Fermi averaged anomalous density Kp(X), EPI V^ S (X) and 
their product, in arbitrary units. 



All the functions are drawn in arbitrary units, for the graphical convenience. It is seen 
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that the anomalous density matrix falls rapidly down outside the slab. This decrease almost 
compensate the surface increase of the EPI in Eq. (|57jl. Indeed, although the product function 
inside the slab is less than at the surface, the difference is not significant. 

To examine a relevance of our previous approximate calculations of A, it is worth to consider 
the so-called gap-shape function which was found in [Tj with the help of a generalization of the 
KKC method j3] to non-homogeneous systems. It was defined as follows: 



X(xi,x 2 ] k\) = A(x 1 ,x 2 ] k 2 ± )/A F (0). 
Obviously, it can be represented as the sum 

x(xx,x 2 ; k]_) = Xi(X) gi(k 2 ±} x), 



(40) 



(41) 



with components 

Xl (X)=A i (X)/A F (0). 
Let us again introduce the Fermi averaged gap-shape function 

XF (X) = 5>(X)^f (X)). 



(42) 



(43) 



This quantity is displayed in Fig. 8 and Fig. 9 at different Eq and K max , correspondingly. 



1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 




Figure 8: The Fermi averaged gap-shape function xf(X) f° r different Eq. 

One can see that the gap-shape function is much less sensitive to a variation of these param- 
eters than the gap itself. In fact, all the sensitivity of the gap is related to the normalization 
factor of Ap(0). Thus, our previous calculations of the gap-shape function with the use of 
K max =60 fm -1 have sufficiently good accuracy. Such high sensitivity of the gap value to tiny 
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Figure 9: The Fermi averaged gap-shape function Xf(X) for different K max . 

variations of the calculation details shows a strong surface- volume interplay. In other words, it 
implies a significant proximity effect. 

The results of the solution of the gap equation for L=8 fm are presented in Fig. 10 (the 
separate components Aj(X)) and in the Table 1 for the diagonal matrix elements A nn (kj_=0), 
n corresponding to the states with the positive parity. 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 



X, fm 

Figure 10: The components of the gap Aj(X). 
These are the matrix elements which could be related to those in heavy nuclei over the 
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single-particle s-states. The typical value of the matrix elements A nn is of the order of 1 
MeV. The only exception is the one corresponding to very small energy e n = — 1.09 MeV. 
The reason is quite trivial. Indeed, the wave function has very long "tail" outside the slab, 
therefore the weight of y n (x) 2 in the integral of A nn is rather small in the region where the 
components Aj(X) are large. The average value of the gap in Table 1 is A=0.85 MeV (and 
= 1.00 MeV, if to exclude the state with the small energy). The self-consistent value of the 
chemical potential /j,— — 7.96 MeV is very close to the one, /x = — 8 MeV, without pairing. The 
difference S/jl— — 0.04 MeV is in agreement with the standard estimate S/j, ~ — A 2 /e F - 

Table 1. Diagonal matrix elements A nn (k ± =0) in the slab with the width parameter L=8 fm. 



e„, MeV 


A nn , MeV 


-49.06 


1.14 


-42.16 


1.11 


-30.22 


0.98 


-15.22 


0.76 


-1.09 


0.26 



Let us now examine dependence of the results on the value of the width parameter L. The 
Fermi averaged gap functions for different L are displayed in Fig. 11. 




X, fm 



Figure 11: The Fermi averaged gap A F (X) for different L. 

The values of L=6.2 fm and L=7.4 fm are chosen to imitate nuclei of the tin and lead 
regions, correspondingly. The average diagonal matrix elements of A(k±=0) are A=1.08 MeV 
for Pb and A=1.06 MeV for Sn. In the latter case, again there is a small energy level with a 
small value of A nn . If it will be excluded from the analysis, one obtains A=1.26 MeV. In any 
case, these values are rather close to those in atomic nuclei which are about 1 MeV for the Pb 
isotopes and about 1.4 MeV, for Sn isotopes. As it is seen in Fig. 11, with increase of L the 
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surface maximum of the function Ap(Jf) becomes lower. At the same time, the central value 
Ap(X=0) becomes smaller. The corresponding values are given in Table 2. At the same time, 
the ratio a 



Ap iax I Af(0) depends on L very smoothly. 



Table 2. Dependence of the cental values Ap(X=0) on the width parameter L. 



L, fm 


A F (A=0) ,MeV 




a 


6.2 


0.45 


1.20 


2.67 


7.4 


0.46 


1.09 


2.37 


8.0 


0.41 


0.94 


2.29 


14.0 


0.36 


0.88 


2.44 


20.0 


0.30 


0.77 


2.57 


oo 


0.22 







The Fermi averaged gap-shape functions at different L are displayed in Fig. 12. In agreement 
with Table 2, the values of the surface maximum are almost L-independent. Fig. 12 is in a 
good agreement with the analogous one in PP obtained with the help of the KKC method. 




10 12 14 16 18 20 22 
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Figure 12: The Fermi averaged gap-shape function xf{X) for different L. 



4 Conclusions 

We solved the gap equation in the 1 S'o-channel for a nuclear slab with the separable 3x3 
representation [21 El of the Paris potential. A two-step method was used in which the gap 
equation is solved in the model space in terms of the EPI which is found in the complementary 
subspace. The model space is chosen sufficiently wide to make it possible to neglect the pairing 
effects in the complementary subspace. In this case, the EPI obeys the Bethe-Goldstone-type 
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equation which was solved with the help of the so-called local-potential approximation which, 
as it was shown in £Q, has very high accuracy. The absolute value of the gap A turned out to 
be very sensitive to the cutoff K max in the momentum space in the equation for the EPI. It is 
necessary to take K max =160 -v- 180 fm _1 to guarantee 1% accuracy for A. 

The gap equation itself is solved directly, without any additional approximations. The 
solution shows the surface enhancement of the gap A which was earlier found ^3] with an 
approximate consideration of the gap equation. A significant surface-volume interplay was 
observed also, thereby a kind of the proximity effect was found. The diagonal matrix elements 
of A which could be corresponded to those in heavy atomic nuclei turned out to be rather close to 
empirical values. Of course, there are various corrections to the BCS approximation. In nuclear 
matter, they came mainly from the screening of the EPI due to particle-hole excitations [15-18] 
and the dispersive self-energy effects [HI]- I n the nuclear slab system under consideration, just 
as in real atomic nuclei, they are partially suppressed due to the enhanced contribution of the 
surface region. Evidently, the most important correction to the BCS approximation is the one 
due to the surface vibrations jB] . In our calculations, the room for all the additional corrections 
to the BCS approximation is of order of 20%. This contradicts to some extent to results of 
Ref . EU HI] where the microscopic gap equation was solved directly for the spherical nucleus 
120 Sn with Argonne NN-force v 14, and only about 50% of the empirical value of A was obtained. 
The additional 50% was explained in [6 with the surface vibration contributions. We suppose 
that the main reason in the difference of our results of solving the BCS gap equation and those 
of Ref. jU] is in value of the effective nucleon mass m* which was chosen in [6 as the density 
dependent one varying from m*=0.7m inside the nucleus to m*=m outside. We used the value 
of m*=m which is, in our opinion, closer to the empirical effective mass in atomic nuclei. As 
to the surface vibration contribution itself, possibly, it is partially overestimated in [S| due to 
neglecting the so-called non-pole diagrams. This point is discussed in more detail in pQ. The 
problem of the consistent account of the non-pole diagrams in the pairing problem, along the 
line of P, is waiting for its solution. 

5 Acknowlegements 

The authors are highly thankful to S.T. Belyaev and S.V. Tolokonnikov for helpful discussions. 
This research was partially supported by the Grant NS-1885.2003.2 of The Russian Ministry 
for Science and Education. 



16 



References 

[I] M. Baldo, U. Lombardo, E. E. Saperstein, M. V. Zverev, Phys. Rep. 391 261 (2004). 
[2] J. Haidenbauer, W. Plessas, Phys. Rev. C 30 (1984) 1822. 

[3] J. Haidenbauer, W. Plessas, Phys. Rev. C 32 (1985) 1424. 

[4] M. Lacombe, B. Loiseaux, J. M. Richard, R. Vinh Mau, J. Cote, D. Pires and R. de Tourreil, 
Phys. Rev. C 21 (1980) 861. 

[5] V. A. Khodel, V. V. Khodel and J. W. Clark, Nucl. Phys. A 598 (1996) 390. 

[6] F. Barranco, R. A. Broglia, G. Colo, G. Gori, E. Vigezzi, and P. F. Bortignon. 
Eur. Phys. J. A 21 57 (2004). 

[7] M. Baldo, J. Cugnon, A. Lejeune, U. Lombardo, Nucl. Phys. A 515 (1990) 409. 

[8] E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. Schaeffer, Nucl. Phys. A 627 (1997) 
710. 

[9] V. A. Khodel and E. E. Saperstein, Phys. Rep. 92 (1982) 183. 

[10] G. E. Brown, J. H. Gunn and P. Gould, Nucl. Phys. 46 (1963) 598. 

[II] A. B. Migdal, Theory of finite Fermi systems and applications to atomic nuclei (Wiley, 
New York, 1967). 

[12] P. Ring, P. Schuck, The nuclear many-body problem (Springer, Berlin, 1980). 

[13] M. Baldo, U. Lombardo, E. E. Saperstein and M. V. Zverev, Nucl. Phys. A 628 (1998) 
503. 

[14] M. Baldo, M. Farine, U. Lombardo, E. E. Saperstein, P. Schuck and M. V. Zverev, 
Eur. Phys. J. A 18 (2003) 17. 

[15] A. Schwenk, B. Friman and G. E. Brown, Nucl. Phys. A 703 (2002) 745. 

[16] T. L. Ainsworth, J. Wambach, and D. Pines, Phys. Lett. B 222 (1989) 173. 

[17] J. Wambach, T. L. Ainsworth, and D. Pines, Nucl. Phys. A 555 (1993) 128. 

[18] H.-J. Schulze, J. Cugnon, A. Lejeune, M. Baldo and U. Lombardo, Phys. Lett. B 375 
(1996) I. 

[19] M. Baldo and A. Grasso, Phys. Lett. B 485 (2000) 115. 

[20] S. G. Kadmensky, and P. A. Luk'yanovich, Sov. J. Nucl. Phys. 49 (1989) 238. 
[21] A. V. Avdeenkov, S. P. Kamerdzhiev, JETP Lett. 69 (1999) 715. 



17 



